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Abstract 



We calculate the third coefficient of the lattice /^-function associated with 
the Wilson formulation for both gauge fields and fermions. This allows us to 
evaluate the three-loop correction (linear in g^) to the relation between the 
lattice A-parameter and the bare coupling go, which is important in order 
to verify asymptotic scaling predictions. Our calculation also leads to the 
two-loop relation between the coupling renormalized in the MS scheme and 
9o- 

The original version of this paper contained a numerical error in one of 
the diagrams, which has now been corrected. The calculations, as well as the 
layout of the paper have remained identical, but there are some important 
changes in the numerical results. 
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I. INTRODUCTION 



In lattice formulations of asymptotically free field theories, such as QCD, the relevant 
region for the continuum physics is that where scaling is verified. Scaling in the continuum 
limit requires that all RG-invariant dimensionless ratios of physical quantities approach 
their continuum value with non-universal corrections which get depressed by integer powers 
of the inverse correlation length. This can be equally stated by introducing a lattice scale 
A L , which is a particular solution of the RG equation 
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(a is the lattice spacing), i.e. 
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and requiring that the ratio of any RG invariant quantity to the appropriate power of Al 
approach a constant in the continuum limit go — > 0. /3l(<7o) is the lattice /3-function 
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which tells us how the bare coupling and the lattice spacing must be changed in order to 
keep physical quantities fixed. 

In the asymptotic region g — > one may perform a weak coupling expansion of the 
/9-function 



Pl(9o) = -b g 3 - hgl - b L 2 gl + 
where, in SU (N) gauge theory with Nf fermion species, 
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Consequently one has the asymptotic relation 
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where 



Q = 
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The verification of the two-loop behavior of A^ is usually called asymptotic scaling. 
Knowledge of the first correction to the two-loop approximation of A^ is important in order 
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to verify the asymptotic prediction ([?]). We recall that Monte Carlo simulations using the 
Wilson action are actually performed at go ~ 1, and so deviations from the two- loop formula 
might not be negligible. b\ has been recently calculated for the pure gauge theory [0,0. In 
this paper we extend this calculation to the full theory including fermions in the Wilson 
formulation. 

The knowledge of b\ can be also used to improve the pertubative relation between the 
coupling g renormalized in the MS (modified minimal subtraction renormalization) scheme 
and the bare lattice coupling g , which is useful in calculations such as studies concerning 
running couplings (see e.g. Refs. Actually, as we shall see, in order to calculate 

&2 , we will first compute the two- loop relation between g and go . This, supplemented with 
knowledge of the three-loop coefficient of the /3-function in the MS scheme, allows us to 
derive b\ . 

In Sec. [n] we describe the method we employed to calculate b\ . We present the results 
of our lattice perturbative calculations, from which we derive the three-loop coefficient b\ 
associated with the Wilson formulation of both gauge fields and fermions. Sec. [Ill] describes 
in some detail the necessary lattice perturbative calculations we performed in the background 
gauge framework. 



II. PERTURBATIVE CALCULATION OF THE /3-FUNCTION. 

The lattice /3-function depends on the lattice formulation considered. We recall that 
only the first two coefficients of its perturbative expansion are universal, i.e. bo and b%. We 
consider the Wilson formulation for both the pure gauge theory and fermions: 

S l = 4 E T r [1 " U^x)) + £ £(4r + '"/.n)'- •/(•'')'••/(•'■) + 

~\ E E ]$A X ) ( r - 7m) U^x^fix + fi)+i) f {x + pi) (r + 7jU ) U^{x) ] i) f {x)] , (9) 

where U^x) is the usual product of link variables U^(x) along the perimeter of a plaquette 
originating at x in the positive fi-u directions; r is the Wilson parameter; / is a flavor index. 

The lattice /3-function is independent of the renormalized fermionic masses. Indeed, 
since it is dimensionless, its mass dependence may come only through the combination 
x = am, where m is the renormalized fermionic mass. Since lattice gauge theories with 
Wilson fermions are renormalizable Jj]], the lattice /3-function must have a well defined zero 
lattice-spacing limit. Thus any dependence on x, and therefore on m, must disappear in 
that limit. In the presence of fermions with different masses, one may also assume a well 
defined chiral limit in order to complete this argument. So in order to calculate /3z(<7o) we 
can restrict ourselves to the case of Nf fermions with zero renormalized mass. 

The relation between the bare lattice coupling go and the renormalized coupling g defined 
in a generic renormalization scheme is 

9o = Z g (go,a/j,)g, (10) 

where \x indicates a renormalization scale. Z g (g Q , a/i) depends on the renormalization scheme. 
The lattice /^-function can be derived from Z g by 
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(3 L (go) = -g a-^\nZ g (g ,aiJ l ) | MiS . (11) 

Thus the evaluation of b\ would require the calculation of Z g (g ,afi) up to three loops on 
the lattice. 

The lattice calculation can be simplified by exploiting the following considerations. Let 
us choose the MS procedure as renormalization scheme in which the renormalized coupling 
g is defined. The corresponding /^-function 

(3(g) = H-Jt |barc quantities= ~b g 3 ~ hg 5 - b 2 g 7 + ••• (12) 

is known to three loops. By writing the MS /3-function in the form 

/3(g) = ln Zg(.9o, a H) \a,g , (13) 



and by comparing it with Eq. (T\\), one can easily derive the relation 

P L (go) = U- 9 2 o^2^K) Z^Z; 1 ), (14) 

which is valid to all orders of perturbation theory. We write Zg as 

Z g (g , an) 2 = 1 + L (afi)g% + L^a/i)^ + 0(gl) (15) 

where L (x) = 2b \iax + l and L\{x) = 2b 1 \nx + l\. The constant Z is related to the 
ratio of the A parameters associated with the particular lattice regularization and the MS 
renormalization scheme: 

/ = 26 ln(A L /A M g). (16) 

Its value as a function of N and Nf can be inferred from Refs. |8|||. It has been obtained 
with improved accuracy in Ref. |]J (pure gluons) and in the present paper (fermionic case). 
For r = 1 

l = - 0.16995599iV + 0.0066960(W f . (17) 
8N 

By expanding Eq. (|Hp in powers of g 2 ., one finds the well-known result that the coefficients 
bo and b\ are the same in (3 L (go) and (3{g). Furthermore 

&£ = &2-&iZo + Wi- (18) 



Thus, since b 2 is known flO 
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the evaluation of b% requires only a two-loop calculation on the lattice, i.e. the calculation 
of the constant l\. 

The computation of the MS renormalization constant Z g is easier in the background field 
gauge. In fact, this renormalization constant has a simple relationship with the background 
field renormalization constant Za fTTH , 



Z A (g ,a/j l )Z g (g ,aiif 



(20) 



Liischer and Weisz |12j have shown that pure lattice gauge theory with a background gauge 
field is renormalizable to all orders in perturbation theory. No additional counterterms are 
required besides those already needed in the absence of a background field. Their argument, 
based on renormalizability of pure lattice gauge theory, BRS, background gauge and shift 
symmetries of the lattice functional integral, can be extended to full lattice QCD in the 
Wilson formulation. An essential point is the renormalizability of lattice gauge theory with 
Wilson fermions proved by Reisz to all orders in perturbation theory [0. 

As a consequence of the relation (|20|) , in order to calculate Z g we need to compute only 
the two-loop self-energy of the background field on the lattice. In the background field 
formulation |l3j the links are written as 



U,(x) = U*{z)U£{z), 



Ufa) 



(21) 



where Qn(x) = T c Q c ^(x) and A^x) = T c y4^(x) are the quantum and background fields 
respectively. The following gauge-fixing term preserves gauge invariance of the background 
field 



S 9f = A oEE Tr ' D / I Qv( X ) D u Qu{x), 
H,v x 

D-Q u (x) = U*~\x - fx)Q u (x - p,)U*(x -p.)- Q v {x). 



(22) 
(23) 



We worked in the Feynman gauge, Ao = 1. This gauge fixing produces the following Fadeev 
Popov action for the ghosts fields u and uJ 



S gh = 2EE Tr {D+oj(x)y(D+u;(x) + ig [Q,(x),u(x)] + \ig Q [Q,(x), D+u(x) 



12 #0 [Ql 



x 



QJx),D+u(x) 



+ ■ ■ 



Dtu(x) = U«(x)cu(x + fim-\x) - uj{x). 



(24) 
(25) 



Finally the change of integration variables from links to vector fields yields a jacobian that 
can be rewritten as a new term S m in the action: 



(26) 



In all the expansions above we have written only the relevant terms for our two-loop com- 
putation. 
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The full action is therefore 

S = Sl + S g f + S g h + S m . (27) 

The vertices needed in perturbation theory are obtained as usual by expanding the expo- 
nential of the action S in powers of the fields. 

As in Ref. [I], we set the following notation. We rewrite the renormalized one-particle 
irreducible two-point functions of the background and quantum fields as 

K\p)t = (<W 2 - PyPu) (1 " "r(p)) l9\ (28) 
r Q R Q (p)t„ = S ab [(W 2 " Pfr) (1 - MP)) + Ap^P,] , (29) 

oo 

MP) = T,9 21 4\P), (30) 
1=1 

oo 

«*(p) = EA4°(p)- (si) 
i=i 



Correspondingly on the lattice 

E rf A {p)t = -S ab 3? [1 - u(p)] /gl (32) 
E ^ Q l Q (pC = ^ ah P 2 1 3 (1 " "(p)) + A ] , (33) 

oo 

Kp) = E^ W (p), (34) 



2=1 

oo 

w(p) = E^ W (p). (35) 
i=i 

The bare and renormalized functions are related by 

[1 - u R (p, fi, g)} = Z A [1 - u(p, a, g )] , (36) 

therefore 

2 _ 1 - v(p, a, g ) 

In order to compare continuum against lattice expressions, one also needs to renormalize 
the gauge parameter according to A = Zq\ |TI| , where Zq is the renormalization constant 
of the quantum field. Zq is needed only to one loop and can be evaluated by imposing 

[1 - u R {p, fi, g)\ = Z Q [1 - u(p, a, g )} . (38) 

One can fix Ao = 1, but then one must take A = Zq as corresponding renormalized gauge 
parameter. 

The MS renormalized functions necessary for the calculation of Z g to two loops are 
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Details of the lattice calculations will be given in the next section. Here we just list the 
results. We report only the fermionic contributions for r = 1 which must be added to the 
pure gauge functions calculated in Refs. piH . 
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ln(aV) + &i/ 
ln(a 2 p 2 ) +hf 



— | 3N - — ) ln(aV) + k 2f — + k 3f N 

(16tt 2 ) 2 V N ) K F J f N 6f 



(42) 
(43) 
• (44) 



k lf = -0.013732194(5), 
k 2f = 0.0011877(14), 
k 3f = -0.0013617(16). 



(45) 



The origin and meaning of the errors in the above results will be explained in the next 
section. 

We have now what we need to calculate Z g to two loops using Eq. ([J7|). The resulting 
two- loop constant 1%, cf. Eq. flT5[) , is given by[] 



h 



128N 2 



+ 0.018127763 - 0.0079101185iV 2 



-0.0011967(14)— + 0.0009998(16)iV 



Through Eq. (jT8f ) one can then obtain b% ■ For example for N = 3 we find 
h\ = -0.0015998323 + 0.0000799(4)^ - 0.00000605(2)]^. 



(46) 



(47) 



Knowledge of b\ allows us to evaluate the coefficient q of the linear correction to asymp- 
totic scaling, cf. Eq. @. For iV = 3, we obtain: 



1 Here and below, results pertaining to the pure gluonic case can be obtained with higher precision 
from Ref. Iffll. 
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q ~ 0.1896 (N f = 0), q ~ 0.2160 (JV/ = 2), g ~ 0.2355 (JV) = 3). (48) 



These numbers show that for g Q ~ 1, which is a typical value of the bare coupling where 
simulations are nowadays performed, the linear correction to asymptotic scaling cannot be 
ignored. 

Finally, the relation between the MS coupling a = g 2 /(4ir) and «o = 5 , o/(^ 7r ) can ^ e 
easily read from Eq. 



a = «o 



(49) 



where, for r 



di(x) = —4rrL (x) 

— (-N - -Nf] \nx - ^ + 2.13573007iV - 0.08414443(8)^/, 



2tt V 3 



(50) 



and 
d 2 {x) = {An) 2 



<ii(x) 2 — 



247T 2 



Lx(x) 



UN 2 - N f ( 13iV - 



A' 



In x 
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- 2.8626216 + 1.2491158A^ + N 



/ 



0.18898(22)— - 0.15789(26)^ 



(51) 



For comparison we have also estimated b\ for r = 0.5 and r = 2.0 . For A^ = 3, we obtain 



&2 — b 2 |jv /= o 



+0.000092N/ - 0.0000033N 



/ 



(52) 



for r = 0.5, and 



^2 ~~ ^2 |iV>=0 



+0.0000079N/ - 0.0000039N 



/ 



(53) 



for r = 2.0 . 



III. THE CALCULATION IN LATTICE PERTURBATION THEORY 

i) Preliminaries 

In this Section we proceed to describe the technical aspects in our calculation of the 
quantities u^(p), u^ 2 \p). The fermionic contributions in cu^(p) are identical to those in 
i/W(p). 

Two diagrams containing fermions contribute to u^(p), shown in Figure 1. There are 
18 two- loop fermionic diagrams contributing to u^(p), as well as two diagrams containing 
an insertion of the one-loop fermion mass counterterm; these are shown in Figure 2. 

The algebra involving lattice quantities was performed using a symbolic manipulation 
package in Mathematica, developed by us in recent years. For the purposes of the present 
work, this package was extended to include fermions. 



8 



A first, relatively brief, part in the evaluation of the diagrams is the contraction. This 
is done completely automatically for diagrams with an arbitrary number of loops, once the 
types of vertices, as well as the "incidence matrix" for the diagram are specified. This step 
also includes: Complete reduction of color structures, Dirac matrices and tensor structures; 
exploiting permutation symmetry and lattice rotational invariance to keep the size of the 
expression down to a minimum; use of trigonometry to arrive at a canonical form. The 
resulting expression is a rational function of sines and cosines of the external (p) and internal 
(q, k) momenta. 

ii) Extracting the external momentum 

The next task is to make explicit the functional dependence of each diagram on p. The 
two-loop amplitude u^(p) can be written as: 

uW(p) = vW(p)\ Nf=0 + 52v i (p) (54) 

i 

(the index i runs over diagrams with fermions shown in Fig. 2), where, generally, 

| / o^p 2 \ 2 In a 2 p 2 j 

ap 2 Vi(p) = co,i + c M a 2 + «V j c 2 ,i j—^ j + c 3|i ^ 2 + c 4>i > + 0((ap) A ) (55) 

(p 2 = 4 ^ sin 2 (p A1 /2)). The dependence of c n ^ of N,Nf is: 

c n>i =[c^ l) /N + c^N]N f (56) 

We note in passing that certain diagrams are infrared convergent only when taken in pairs: 
(7, 11), (8, 18), (9, 17); we have evaluated these accordingly, with due care taken to avoid 
divergences in intermediate results. 

To extract the p-dependence, we first isolate the superficially divergent terms; these are 
responsible for the double logarithms. There are relatively few such terms, and in the pure 
gluonic case they all have been tabulated 0. We can use these tables also in diagrams with 
fermions, applying successive subtractions of the type: 

q 2 q 2 ~^ (g 2 q 2 ) ^ ^ 

where if is the inverse fermionic propagator: q 2 = (q 2 r/2) 2 + S/i sm2 Qyi • This leads to the 
tabulated expressions plus a series of superficially convergent terms. 

The rest of the expression may still contain sub divergences, which can give rise to single 
logarithms. These are handled using successive subtractions of the type: 



1 / 1 



-,2 T 2 1 —2 

q + k k \q + k 



(58) 



The first term on the right hand side then yields factorized one-loop expressions, whose 
p-dependence is easily extracted, while the second term is more convergent. Finally, all 
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remaining terms, containing no (sub-) divergences can be evaluated by Taylor expansion in 
ap. 

In our code, judicious choices for the right set of subtractions are applied automatically to 
each term. While the extraction of divergences is not particularly complicated conceptually, 
it leads to a great proliferation in the resulting number of terms. As an example, diagram 16 
by itself leads to ~100 types of expressions; each type must then be numerically evaluated on 
its own (because it contains a different set of subtracted propagators) and contains typically 
up to some hundreds of terms. 

Hi) Numerical integration 

At this stage, all the above types of expressions no longer contain p and must be numer- 
ically integrated over the internal momenta. The integration is done in momentum space 
over finite lattices; an extrapolation to infinite size is then performed (see below). 

For expressions with only gluonic propagators a coordinate space method was proposed 
li| , yielding high precision results. Recently, this method was also used for some study 



in 

cases involving fermions |I5[]. For the present calculation, involving fermionic and bosonic 
propagators with multiple subtractions, the coordinate space method becomes a bit more 
complicated; we may return to it in the future, if the need arises for more precision than 
presented in the present paper. 

Fortran code for the numerical evaluation of each type of expression is then created 
automatically by our programs. The code is highly optimized: It avoids redundant evalua- 
tion of common subexpressions, of symmetric regions in momentum space, of r-independent 
subexpressions, etc. The code can be run for different values of r and lattice size L . For 
the present paper we used L < 30; we present r = 1 and, for comparison, we also estimated 
r = 0.5, r = 2.0. 

Extrapolation to infinite lattice size is of course a source of systematic error. To esti- 
mate this error, our procedure carries out automatically the following steps: First, different 
extrapolations are performed using a broad spectrum of functional forms of the type: 

^e^L^ (In LY (59) 

For the k th such extrapolation, a deviation dk is calculated using several criteria for quality of 
fit. Finally, these deviations are used to assign weights d^ 2 j (J2k d^ 2 ) to each extrapolation, 
producing a final value together with the error estimate. We have checked the validity of 
these estimates in cases where the exact answer was known (see also below), finding the 
estimate to be always correct. 

iv) Cancellations and cross checks 

Several constraints exist on the coefficients c n) j. We have used these constraints as 
verifications both on the algebraic expressions and on the numerical results: 

co,j: Gauge invariance requires 

E C <M = 0- (60) 

i 

We checked this property in three ways: Firstly, a formal, Ward-identity type derivation 
was performed, in which vertices with background fields at zero momentum were written in 
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terms of appropriate derivatives of inverse propagators. From this we find not only Eq. (|60| ) . 
but also some additional constraints: 



2 Co,l + Co,3 — 

Co,l + C ,3 + C ,6 + C ,8 = 

(C ,4 + C , 5 )|7V2=2 = 

Co,7 + C ,9 + Co,H + |c ,l6 + C ,17 = 

Co,13 = 

c 0,19 + c 0,20 — 

(c 0l i2 + c ,i 8 )(iV 2 - 2) + c ,i 4 (iV 2 - 1) - c ,i 5 (iV 2 - l)(iV 2 - 2) = (61) 

The above are true for any value of the fermion mass m and Wilson parameter r. Secondly, 
all the above identities were checked by algebraic manipulation of the full expression for 
each of the coefficients. Thirdly, we substituted in each identity the numerical results for 
the coefficients, finding in all cases zero within the error estimates. 

c\ i\ The sum of these terms must vanish if Lorentz invariance is to be recovered in the 
continuum limit. Only diagrams 9 and 17 give nonvanishing contributions, and we checked 
in the same ways as above that their sum is zero. 

C2,f These coefficients must coincide with those of the continuum. We checked that this is 
so: 

11 4 5 liV 2 -l 

02,15 = g-j^V) , c 2 , 16 = -NN f , c 2)17 = --NN f , c 2 , 18 = 3— j^V> • (62) 

For all other diagrams: C2,i = 0. 

c 3i : The total contribution for single logarithms must coincide with the continuum result: 

S^is^-^w (63) 

i 

Again, this was checked both algebraically and numerically. 

The results for each diagram are presented in Tables 1 and 2. Diagrams not appearing 
in these Tables give vanishing contributions.^] 

For the one- loop amplitude ap u^(p), the individual contributions from the correspond- 
ing two diagrams in Figure 1 are, respectively: 

0.040848920(4) - a 2 p 2 0.013732194(5) + a 2 p 2 ln(a 2 p 2 ) , -0.040848919(5) (64) 

24tt 

In concluding this section, we would like to point out that the procedure outlined here 
applies unchanged to several other interesting cases, including: More complicated actions, 



Compared to the original version, only one entry is different, corresponding to diagram 19. 
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matrix elements of different operators. We hope to address these issues in a future publica- 
tion. 
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TABLES 



TABLE I. Coefficients c ( /> , c\/\ c\ yi ' . r = 1. 



i 


C 0,i 


c 3,i 




1 


-0.00158221542(13) 








3 


0.0031644309(6) 


0.051644463410 


-0.0010637877(8) 


4 


0.00039273(5) 








5 


-0.0005077(9) 





0.00010879(5) 


6 


-0.0086173230(12) 








7+11 








0.00069292(12) 


8 


0.0070351064(17) 


0.258222317052 


-0.0031633707(22) 


12 


0.00069739(6) 








14 


0.00022988(14) 


0.00528244566(5) 


0.0000757(12) 


15 


-0.00060983(5) 


-0.019554587(16) 


0.0001732(6) 


18 


-0.0002019(4) 


0.023786891(14) 


-0.00007910(15) 


19 


-0.010721166(2) 


-0.325714118(9) 


0.004443346(3) 


20 


0.0107211662(6) 









TABLE II. Coefficients cjft, 4 /, c£< . r = 1. 



i 


c (1) 


c (1) 


c {1) 


1 


0.00158221542(13) 








3 


-0.0031644309(6) 


-0.051644463410 


0.0010637877(8) 


4 


-0.00039273(5) 








5 


0.0004503(6) 





-0.00008893(4) 


6 


0.0086173230(12) 








7+11 


0.0014560(8) 





-0.0007427(4) 


8 


-0.0070351064(17) 


-0.258222317052 


0.0031633707(22) 


9+17 


-0.0020085(5) 


0.08132606(7) 


-0.0011335(13) 


12 


-0.00069739(6) 








13 








0.00006908(5) 


14 


-0.00011494(7) 


-0.00264122283(3) 


-0.0000378(6) 


16 


0.0011052(4) 


-0.0517470(17) 


0.0007092(7) 


18 


0.0002019(4) 


-0.023786891(14) 


0.00007910(15) 


19 


0.010721166(2) 


0.325714118(9) 


-0.004443346(3) 


20 


-0.0107211662(6) 
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FIGURES 



FIG. 1. Diagrams with fermionic lines contributing to the one-loop function v^\p). Dashed 
lines ending on a cross represent background gluons. Solid lines represent fermions. 

FIG. 2. Diagrams with fermionic lines contributing to the two-loop function v^ 2 \p). Dashed 
lines represent gluonic fields; those ending on a cross stand for background gluons. Solid lines 
represent fermions. The filled circle is a one-loop fermion mass counterterm. 
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